Increased power by harmonizing structural MRI site differences with the ComBat batch adjustment method in ENIGMA

A common limitation of neuroimaging studies is their small sample sizes. To overcome this hurdle, the Enhancing Neuro Imaging Genetics through Meta-Analysis (ENIGMA) Consortium combines neuroimaging data from many institutions worldwide. However, this introduces heterogeneity due to different scanning devices and sequences. ENIGMA projects commonly address this heterogeneity with random-effects meta-analysis or mixed-effects mega-analysis. Here we tested whether the batch adjustment method, ComBat, can further reduce site-related heterogeneity and thus increase statistical power. We conducted random-effects meta-analyses, mixed-effects mega-analyses and ComBat mega-analyses to compare cortical thickness, surface area and subcortical volumes between 2897 individuals with a diagnosis of schizophrenia and 3141 healthy controls from 33 sites. Specifically, we compared the imaging data between individuals with schizophrenia and healthy controls, covarying for age and sex. The use of ComBat substantially increased the statistical significance of the findings as compared to random-effects meta-analyses. The findings were more similar when comparing ComBat with mixed-effects mega-analysis, although ComBat still slightly increased the statistical significance. ComBat also showed increased statistical power when we repeated the analyses with fewer sites. Results were nearly identical when we applied the ComBat harmonization separately for cortical thickness, cortical surface area and subcortical volumes. Therefore, we recommend applying the ComBat function to attenuate potential effects of site in ENIGMA projects and other multi-site structural imaging work. We provide easy-to-use functions in R that work even if imaging data are partially missing in some brain regions, and they can be trained with one data set and then applied to another (a requirement for some analyses such as machine learning)

Large-scale analysis of structural brain asymmetries in schizophrenia via the ENIGMA consortium

BACKGROUND Left-right asymmetry is an important organizing feature of the healthy brain that may be altered in schizophrenia, but most studies have used relatively small samples and heterogeneous approaches, resulting in equivocal findings. We carried out the largest case-control study of structural brain asymmetries in schizophrenia (N = 11,095), using a single image analysis protocol. METHODS We included T1-weighted data from 46 datasets (5,080 affected individuals and 6,015 controls) from the ENIGMA Consortium. Asymmetry indexes were calculated for global and regional cortical thickness, surface area, and subcortical volume measures. Differences of asymmetry were calculated between affected individuals and controls per dataset, and effect sizes were meta-analyzed across datasets. Analyses were also performed with respect to the use of antipsychotic medication and other clinical variables, as well as age and sex. Case-control differences in a multivariate context were assessed in a subset of the data (N = 2,029). RESULTS Small average differences between cases and controls were observed for asymmetries in cortical thickness, specifically of the rostral anterior cingulate (d = −0.08, pFDR = 0.047) and the middle temporal gyrus (d = −0.07, pFDR = 0.048), both driven primarily by thinner cortices in the left hemisphere in schizophrenia. These asymmetries were not significantly associated with the use of antipsychotic medication or other clinical variables. Older individuals with schizophrenia showed a stronger average leftward asymmetry of pallidum volume than older controls (d = 0.08, pFDR = 9.0 × 10−3). The multivariate analysis revealed that 7% of the variance across all structural asymmetries was explained by case-control status (F = 1.87, p = 1.25 × 10−5). CONCLUSIONS Altered trajectories of asymmetrical brain development and/or lifespan asymmetry may contribute to schizophrenia pathophysiology. Small case-control differences of brain macro-structural asymmetry may manifest due to more substantial differences at the molecular, cytoarchitectonic or circuit levels, with functional relevance for lateralized cognitive processes

What we learn about bipolar disorder from large-scale neuroimaging: Findings and future directions from theENIGMABipolar Disorder Working Group

MRI‐derived brain measures offer a link between genes, the environment and behavior and have been widely studied in bipolar disorder (BD). However, many neuroimaging studies of BD have been underpowered, leading to varied results and uncertainty regarding effects. The Enhancing Neuro Imaging Genetics through Meta‐Analysis (ENIGMA) Bipolar Disorder Working Group was formed in 2012 to empower discoveries, generate consensus findings and inform future hypothesis‐driven studies of BD. Through this effort, over 150 researchers from 20 countries and 55 institutions pool data and resources to produce the largest neuroimaging studies of BD ever conducted. The ENIGMA Bipolar Disorder Working Group applies standardized processing and analysis techniques to empower large‐scale meta‐ and mega‐analyses of multimodal brain MRI and improve the replicability of studies relating brain variation to clinical and genetic data. Initial BD Working Group studies reveal widespread patterns of lower cortical thickness, subcortical volume and disrupted white matter integrity associated with BD. Findings also include mapping brain alterations of common medications like lithium, symptom patterns and clinical risk profiles and have provided further insights into the pathophysiological mechanisms of BD. Here we discuss key findings from the BD working group, its ongoing projects and future directions for large‐scale, collaborative studies of mental illness

Do Rapoport's Rule, Mid-Domain Effect or Environmental Factors Predict Latitudinal Range Size Patterns of Terrestrial Mammals in China?

BACKGROUND: Explaining species range size pattern is a central issue in biogeography and macroecology. Although several hypotheses have been proposed, the causes and processes underlying range size patterns are still not clearly understood. In this study, we documented the latitudinal mean range size patterns of terrestrial mammals in China, and evaluated whether that pattern conformed to the predictions of the Rapoport's rule several analytical methods. We also assessed the influence of the mid-domain effect (MDE) and environmental factors on the documented range size gradient. METHODOLOGY/PRINCIPAL FINDINGS: Distributions of 515 terrestrial mammals and data on nine environmental variables were compiled. We calculated mean range size of the species in each 5° latitudinal band, and created a range size map on a 100 km×100 km quadrat system. We evaluated Rapoport's rule according to Steven's, mid-point, Pagel's and cross-species methods. The effect of the MDE was tested based on a Monte Carlo simulation and linear regression. We used stepwise generalized linear models and correlation analyses to detect the impacts of mean climate condition, climate variability, ambient energy and topography on range size. The results of the Steven's, Pagel's and cross-species methods supported Rapoport's rule, whereas the mid-point method resulted in a hump-shaped pattern. Our range size map showed that larger mean latitudinal extents emerged in the mid-latitudes. We found that the MDE explained 80.2% of the range size variation, whereas, environmental factors accounted for <30% of that variation. CONCLUSIONS/SIGNIFICANCE: Latitudinal range size pattern of terrestrial mammals in China supported Rapoport's rule, though the extent of that support was strongly influenced by methodology. The critical factor underlying the observed gradient was the MDE, and the effects of climate, energy and topography were limited. The mean climate condition hypothesis, climate variability hypothesis, ambient energy hypotheses and topographical heterogeneity hypotheses were not supported

Diagnosis of bipolar disorders and body mass index predict clustering based on similarities in cortical thickness-ENIGMA study in 2436 individuals

AIMS: Rates of obesity have reached epidemic proportions, especially among people with psychiatric disorders. While the effects of obesity on the brain are of major interest in medicine, they remain markedly under-researched in psychiatry. METHODS: We obtained body mass index (BMI) and magnetic resonance imaging-derived regional cortical thickness, surface area from 836 bipolar disorders (BD) and 1600 control individuals from 14 sites within the ENIGMA-BD Working Group. We identified regionally specific profiles of cortical thickness using K-means clustering and studied clinical characteristics associated with individual cortical profiles. RESULTS: We detected two clusters based on similarities among participants in cortical thickness. The lower thickness cluster (46.8% of the sample) showed thinner cortex, especially in the frontal and temporal lobes and was associated with diagnosis of BD, higher BMI, and older age. BD individuals in the low thickness cluster were more likely to have the diagnosis of bipolar disorder I and less likely to be treated with lithium. In contrast, clustering based on similarities in the cortical surface area was unrelated to BD or BMI and only tracked age and sex. CONCLUSIONS: We provide evidence that both BD and obesity are associated with similar alterations in cortical thickness, but not surface area. The fact that obesity increased the chance of having low cortical thickness could explain differences in cortical measures among people with BD. The thinner cortex in individuals with higher BMI, which was additive and similar to the BD-associated alterations, may suggest that treating obesity could lower the extent of cortical thinning in BD

Immunoglobulin, glucocorticoid, or combination therapy for multisystem inflammatory syndrome in children: a propensity-weighted cohort study

Background: Multisystem inflammatory syndrome in children (MIS-C), a hyperinflammatory condition associated with SARS-CoV-2 infection, has emerged as a serious illness in children worldwide. Immunoglobulin or glucocorticoids, or both, are currently recommended treatments. Methods: The Best Available Treatment Study evaluated immunomodulatory treatments for MIS-C in an international observational cohort. Analysis of the first 614 patients was previously reported. In this propensity-weighted cohort study, clinical and outcome data from children with suspected or proven MIS-C were collected onto a web-based Research Electronic Data Capture database. After excluding neonates and incomplete or duplicate records, inverse probability weighting was used to compare primary treatments with intravenous immunoglobulin, intravenous immunoglobulin plus glucocorticoids, or glucocorticoids alone, using intravenous immunoglobulin as the reference treatment. Primary outcomes were a composite of inotropic or ventilator support from the second day after treatment initiation, or death, and time to improvement on an ordinal clinical severity scale. Secondary outcomes included treatment escalation, clinical deterioration, fever, and coronary artery aneurysm occurrence and resolution. This study is registered with the ISRCTN registry, ISRCTN69546370. Findings: We enrolled 2101 children (aged 0 months to 19 years) with clinically diagnosed MIS-C from 39 countries between June 14, 2020, and April 25, 2022, and, following exclusions, 2009 patients were included for analysis (median age 8·0 years [IQR 4·2–11·4], 1191 [59·3%] male and 818 [40·7%] female, and 825 [41·1%] White). 680 (33·8%) patients received primary treatment with intravenous immunoglobulin, 698 (34·7%) with intravenous immunoglobulin plus glucocorticoids, 487 (24·2%) with glucocorticoids alone; 59 (2·9%) patients received other combinations, including biologicals, and 85 (4·2%) patients received no immunomodulators. There were no significant differences between treatments for primary outcomes for the 1586 patients with complete baseline and outcome data that were considered for primary analysis. Adjusted odds ratios for ventilation, inotropic support, or death were 1·09 (95% CI 0·75–1·58; corrected p value=1·00) for intravenous immunoglobulin plus glucocorticoids and 0·93 (0·58–1·47; corrected p value=1·00) for glucocorticoids alone, versus intravenous immunoglobulin alone. Adjusted average hazard ratios for time to improvement were 1·04 (95% CI 0·91–1·20; corrected p value=1·00) for intravenous immunoglobulin plus glucocorticoids, and 0·84 (0·70–1·00; corrected p value=0·22) for glucocorticoids alone, versus intravenous immunoglobulin alone. Treatment escalation was less frequent for intravenous immunoglobulin plus glucocorticoids (OR 0·15 [95% CI 0·11–0·20]; p<0·0001) and glucocorticoids alone (0·68 [0·50–0·93]; p=0·014) versus intravenous immunoglobulin alone. Persistent fever (from day 2 onward) was less common with intravenous immunoglobulin plus glucocorticoids compared with either intravenous immunoglobulin alone (OR 0·50 [95% CI 0·38–0·67]; p<0·0001) or glucocorticoids alone (0·63 [0·45–0·88]; p=0·0058). Coronary artery aneurysm occurrence and resolution did not differ significantly between treatment groups. Interpretation: Recovery rates, including occurrence and resolution of coronary artery aneurysms, were similar for primary treatment with intravenous immunoglobulin when compared to glucocorticoids or intravenous immunoglobulin plus glucocorticoids. Initial treatment with glucocorticoids appears to be a safe alternative to immunoglobulin or combined therapy, and might be advantageous in view of the cost and limited availability of intravenous immunoglobulin in many countries. Funding: Imperial College London, the European Union's Horizon 2020, Wellcome Trust, the Medical Research Foundation, UK National Institute for Health and Care Research, and National Institutes of Health

Widespread white matter microstructural differences in schizophrenia across 4322 individuals:Results from the ENIGMA Schizophrenia DTI Working Group

The regional distribution of white matter (WM) abnormalities in schizophrenia remains poorly understood, and reported disease effects on the brain vary widely between studies. In an effort to identify commonalities across studies, we perform what we believe is the first ever large-scale coordinated study of WM microstructural differences in schizophrenia. Our analysis consisted of 2359 healthy controls and 1963 schizophrenia patients from 29 independent international studies; we harmonized the processing and statistical analyses of diffusion tensor imaging (DTI) data across sites and meta-analyzed effects across studies. Significant reductions in fractional anisotropy (FA) in schizophrenia patients were widespread, and detected in 20 of 25 regions of interest within a WM skeleton representing all major WM fasciculi. Effect sizes varied by region, peaking at (d=0.42) for the entire WM skeleton, driven more by peripheral areas as opposed to the core WM where regions of interest were defined. The anterior corona radiata (d=0.40) and corpus callosum (d=0.39), specifically its body (d=0.39) and genu (d=0.37), showed greatest effects. Significant decreases, to lesser degrees, were observed in almost all regions analyzed. Larger effect sizes were observed for FA than diffusivity measures; significantly higher mean and radial diffusivity was observed for schizophrenia patients compared with controls. No significant effects of age at onset of schizophrenia or medication dosage were detected. As the largest coordinated analysis of WM differences in a psychiatric disorder to date, the present study provides a robust profile of widespread WM abnormalities in schizophrenia patients worldwide. Interactive three-dimensional visualization of the results is available at www.enigma-viewer.org.Molecular Psychiatry advance online publication, 17 October 2017; doi:10.1038/mp.2017.170

Call for emergency action to restore dietary diversity and protect global food systems in times of COVID-19 and beyond: Results from a cross-sectional study in 38 countries

Background: The COVID-19 pandemic has revealed the fragility of the global food system, sending shockwaves across countries\u27 societies and economy. This has presented formidable challenges to sustaining a healthy and resilient lifestyle. The objective of this study is to examine the food consumption patterns and assess diet diversity indicators, primarily focusing on the food consumption score (FCS), among households in 38 countries both before and during the first wave of the COVID-19 pandemic. Methods: A cross-sectional study with 37 207 participants (mean age: 36.70 ± 14.79, with 77 % women) was conducted in 38 countries through an online survey administered between April and June 2020. The study utilized a pre-tested food frequency questionnaire to explore food consumption patterns both before and during the COVID-19 periods. Additionally, the study computed Food Consumption Score (FCS) as a proxy indicator for assessing the dietary diversity of households. Findings: This quantification of global, regional and national dietary diversity across 38 countries showed an increment in the consumption of all food groups but a drop in the intake of vegetables and in the dietary diversity. The household\u27s food consumption scores indicating dietary diversity varied across regions. It decreased in the Middle East and North Africa (MENA) countries, including Lebanon (p \u3c 0.001) and increased in the Gulf Cooperation Council countries including Bahrain (p = 0.003), Egypt (p \u3c 0.001) and United Arab Emirates (p = 0.013). A decline in the household\u27s dietary diversity was observed in Australia (p \u3c 0.001), in South Africa including Uganda (p \u3c 0.001), in Europe including Belgium (p \u3c 0.001), Denmark (p = 0.002), Finland (p \u3c 0.001) and Netherland (p = 0.027) and in South America including Ecuador (p \u3c 0.001), Brazil (p \u3c 0.001), Mexico (p \u3c 0.0001) and Peru (p \u3c 0.001). Middle and older ages [OR = 1.2; 95 % CI = [1.125–1.426] [OR = 2.5; 95 % CI = [1.951–3.064], being a woman [OR = 1.2; 95 % CI = [1.117–1.367], having a high education (p \u3c 0.001), and showing amelioration in food-related behaviors [OR = 1.4; 95 % CI = [1.292–1.709] were all linked to having a higher dietary diversity. Conclusion: The minor to moderate changes in food consumption patterns observed across the 38 countries within relatively short time frames could become lasting, leading to a significant and prolonged reduction in dietary diversity, as demonstrated by our findings

In vivo hippocampal subfield volumes in bipolar disorder—A mega-analysis from The Enhancing Neuro Imaging Genetics through Meta-Analysis Bipolar Disorder Working Group

The hippocampus consists of anatomically and functionally distinct subfields that may be differentially involved in the pathophysiology of bipolar disorder (BD). Here we, the Enhancing NeuroImaging Genetics through Meta‐Analysis Bipolar Disorder workinggroup, study hippocampal subfield volumetry in BD. T1‐weighted magnetic resonance imaging scans from 4,698 individuals (BD = 1,472, healthy controls [HC] = 3,226) from 23 sites worldwide were processed with FreeSurfer. We used linear mixed‐effects models and mega‐analysis to investigate differences in hippocampal subfield volumes between BD and HC, followed by analyses of clinical characteristics and medication use. BD showed significantly smaller volumes of the whole hippocampus (Cohen's d = −0.20), cornu ammonis (CA)1 (d = −0.18), CA2/3 (d = −0.11), CA4 (d = −0.19), molecular layer (d = −0.21), granule cell layer of dentate gyrus (d = −0.21), hippocampal tail (d = −0.10), subiculum (d = −0.15), presubiculum (d = −0.18), and hippocampal amygdala transition area (d = −0.17) compared to HC. Lithium users did not show volume differences compared to HC, while non‐users did. Antipsychotics or antiepileptic use was associated with smaller volumes. In this largest study of hippocampal subfields in BD to date, we show widespread reductions in nine of 12 subfields studied. The associations were modulated by medication use and specifically the lack of differences between lithium users and HC supports a possible protective role of lithium in BD

Using structural MRI to identify bipolar disorders - 13 site machine learning study in 3020 individuals from the ENIGMA Bipolar Disorders Working Group.

Bipolar disorders (BDs) are among the leading causes of morbidity and disability. Objective biological markers, such as those based on brain imaging, could aid in clinical management of BD. Machine learning (ML) brings neuroimaging analyses to individual subject level and may potentially allow for their diagnostic use. However, fair and optimal application of ML requires large, multi-site datasets. We applied ML (support vector machines) to MRI data (regional cortical thickness, surface area, subcortical volumes) from 853 BD and 2167 control participants from 13 cohorts in the ENIGMA consortium. We attempted to differentiate BD from control participants, investigated different data handling strategies and studied the neuroimaging/clinical features most important for classification. Individual site accuracies ranged from 45.23% to 81.07%. Aggregate subject-level analyses yielded the highest accuracy (65.23%, 95% CI = 63.47-67.00, ROC-AUC = 71.49%, 95% CI = 69.39-73.59), followed by leave-one-site-out cross-validation (accuracy = 58.67%, 95% CI = 56.70-60.63). Meta-analysis of individual site accuracies did not provide above chance results. There was substantial agreement between the regions that contributed to identification of BD participants in the best performing site and in the aggregate dataset (Cohen's Kappa = 0.83, 95% CI = 0.829-0.831). Treatment with anticonvulsants and age were associated with greater odds of correct classification. Although short of the 80% clinically relevant accuracy threshold, the results are promising and provide a fair and realistic estimate of classification performance, which can be achieved in a large, ecologically valid, multi-site sample of BD participants based on regional neurostructural measures. Furthermore, the significant classification in different samples was based on plausible and similar neuroanatomical features. Future multi-site studies should move towards sharing of raw/voxelwise neuroimaging data